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Abstract: We show that the complete planar S-matrix of M = 4 super Yang- Mills - 
including all N fc MHV partial amplitudes to all loops - is equivalent to the correlation 
function of a supersymmetric Wilson loop in twistor space. Remarkably, the entire 
classical S-matrix arises from evaluating the correlation function in the self-dual 
sector, while the expansion of the correlation function in powers of the Yang-Mills 
coupling constant provides the loop expansion of the amplitudes. We support our 
proposal with explicit computations of the n particle NMHV and N 2 MHV trees, the 
integrands of the 1-loop MHV and NMHV amplitudes, and the n particle 2-loop 
MHV amplitude. These calculations are performed using the twistor action in ax- 
ial gauge. In this gauge, the Feynman diagrams of the correlation function are the 
planar duals of the usual MHV diagrams for the scattering amplitude. The results 
are presented in the form of a sum of products of dual superconformal invariants 
in (momentum) twistor space, and agree with the expressions derived in the com- 
panion paper [1] directly from the MHV rules. The twistor space Wilson loop is a 
natural supersymmetric generalization of the standard Wilson loop used to compute 
MHV amplitudes. We show how the Penrose- Ward transform can be used to deter- 
mine a corresponding supersymmetrization on space-time and give the corresponding 
superconnection in the abelian case. 
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1. Introduction 

The main purpose of this paper is to show that all n-particle planar amplitudes in 
M = 4 super Yang-Mills are captured by the correlation function of a supersymmetric 
Wilson loop 

<W[C„]> (1.1) 

in (momentum) twistor space. This includes both trees and all loop corrections for 
all N fc MHV partial amplitudes, up to an overall factor of the MHV tree. 

We defer the full explanation of this operator to section 2, but in brief it is a 
supersymmetric generalization of the translation into twistor space of a space-time 
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Figure 1: The nodal curve C n in twistor space that corresponds to the null polygon in 
space-time. The twistors Z^ represent the edges of the space-time polygon, while the lines 
(Zi, Zi-x) correspond to the space-time vertices Xi. 

Wilson loop associated to a null polygon. The Wilson loop W[C n ] is a natural holo- 
morphic analogue of the standard space-time definition of the trace of the holonomy 
of a connection around a loop: it computes the trace of the holonomy of the process 
of finding holomorphic frames of an almost complex bundle around a nodal curve 
C n in supertwistor space CP 3 ' 4 . This curve is a polygon whose 'edges' are Riemann 
spheres (complex projective lines) dual to the vertices of an n-sided null polygon 
C' n in dual conformal space-time (see figure 1). For much of the discussion here, 
all that is needed of the geometry of the twistor correspondence is that points in 
space-time correspond to complex projective lines (CP 1 s) in twistor space, in such a 
way that intersection of lines in twistor space is equivalent to the points lying on a 
null geodesic. See e.g. [2, 3, 4] for introductions to twistor geometry. 

The relation to scattering amplitudes comes from treating the twistor space to 
be the space of the fundamental representation of the dual superconformal group [5] , 
rather than the usual superconformal group. This dual group acts on the space of 
'region momenta' of a scattering amplitude - defined up to translation by Xi — Xj+i = 
Pi where pi is the momentum of the i th particle in the scattering process - in exactly 
the same way as the usual superconformal group acts on space-time. The twistor 
space associated to region momentum space is called momentum twistor space; it 
has the same relation to dual space-time as ordinary twistor space has to usual 
space-time. (See e.g. [6, 7] for more detailed introductions to momentum twistors.) 

Our proposal is motivated by the corresponding proposal in dual space-time that 
null polygonal Wilson loops yield the planar MHV amplitude to all loop orders [8, 
9, 10, 11, 12, 13, 14, 15, 16]. This was in turn motivated by the proposal that taking 
the polygon to bound a minimal surface in AdSs, the area of this surface gives the 
strong coupling limit of the same planar MHV amplitude [17, 18, 19, 20]. However, 
on dual space-time it has not yet been possible obtain non MHV amplitudes in this 
way. 

To support our main conjecture, we use the correlation function (1.1) to obtain 
the NMHV and N 2 MHV tree amplitudes, and the integrands of the MHV and NMHV 
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1-loop and MHV 2-loop amplitudes, each for an arbitrary number of particles. These 
calculations turn out to be straightforward using the twistor action for Af = 4 SYM 
that was introduced in [21, 22, 23]. The twistor action is reviewed in section 3. It 
consists of a holomorphic Chern-Simons theory, corresponding to self-dual M = 4 
SYM, together with an infinite series of MHV vertices that are supported on lines in 
twistor space. Since it lives on a six (real) dimensional space, the twistor action has 
a larger amount of gauge freedom than the usual space-time action, providing access 
to gauges not available on space-time. In particular, in an axial gauge, its Feynman 
diagrams are the MHV diagrams of [24], here being used to compute a correlation 
function rather than a scattering amplitude. 

In the companion paper [1], we translated the MHV formalism for scattering 
amplitudes into momentum twistor space. There, we found that the vertices of an 
MHV diagram are represented simply by 1 on momentum twistor space, while each 
propagator corresponds to an R invariant [5], constructed from the pair of twistors 
associated to each of the two regions separated by the propagator, and an additional 
reference twistor. Remarkably, the diagrams generated by the correlator (1.1) are 
precisely the planar dual graphs of the MHV diagrams for the corresponding planar 
amplitude! In particular, the number V of MHV vertices in the Feynman diagrams 
of the correlation function (1.1) is simply the loop order I of the scattering process. 
Consequently, the entire classical S-matrix of M = 4 SYM comes from the correlation 
function of the supersymmetric Wilson loop operator purely in the holomorphic 
Chern-Simons theory, as we will see in section 4. 

In section 5 we consider corrections to the correlator (1.1) in powers of the Yang- 
Mills coupling. The fact that the diagrams of the twistor Wilson loop (1.1) are dual 
to the usual MHV diagram, with the expression agreeing identically even before 
performing the loop integrals, makes it obvious that the Wilson loop corresponds 
to the scattering amplitude, at least at a computational level. This is in marked 
contrast to the space-time calculation, where the relation between the Wilson loop 
and scattering amplitudes is mysterious, with seemingly very different constructions 
providing identical answers at the final stage (compare [10] and [25], for example). 

Using the twistor action naturally packages the result of the correlator as the 
same sum of products of R invariants as found in [1] from the MHV formalism for 
scattering amplitudes. To compute the loop amplitudes, one must integrate these R 
invariants over copies of the twistor space, two for each loop. Loop integrals over 
twistor space are often tractable, and have been explicitly carried out in various 1- 
and 2-loop cases [26, 27, 28, 29] using the Coulomb branch regulator introduced in [30] 
to treat the infra-red divergences. However, the integrand of the loop amplitude - 
a concept that makes sense only in the planar limit where one may use the cyclic 
ordering to define a loop momentum common to different diagrams - is an interesting 
object in its own right, and has recently been proved [31] to possess the full Yangian 
invariance of planar J\f = 4 SYM, first discovered in scattering amplitudes in [32]. 
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It is therefore intriguing that the correspondence between scattering amplitudes and 
the Wilson loop appears to hold even at the level of the integrand. 

Finally, our twistor space Wilson loop is a supersymmetric formulation of a tran- 
scription of the standard space-time Wilson loop. Although the twistor formulation 
is especially elegant, the results of this paper strongly suggest that it is also possible 
to find a generalisation of the Wilson loop to all N fc MHV amplitudes purely on space- 
time. In section 6 we explain how an off-shell version of the Penrose- Ward transform 
can be used to deduce the superconnection explicitly when it does correspond to one. 
The trace of the holonomy of this superconnection should give rise to an equivalent 
space-time supersymmetric Wilson loop. In space-time, the twistor action most nat- 
urally corresponds to the Chalmers- Siegel action for M = 4 SYM [33, 34]. From 
this point of view, the MHV vertices that give rise to loop amplitudes essentially 
correspond to the Lagrangian insertions considered in [35, 36, 37]. 

The ability to compute correlation functions using the twistor action is not at 
all limited to the specific operator in (1.1). The ideas of this paper will be pre- 
sented in greater generality, with fuller explanations and further elaboration, in the 
forthcoming paper [38]. 

2. Twistor space and Wilson loops 

We begin by defining the operator in (1.1), explaining in general terms why this 
operator is expected to be related to a space-time Wilson loop. 

2.1 The off-shell Penrose- Ward correspondence 

To describe gauge theories on twistor space, we introduce a rank N complex bundle 
E over twistor space that has vanishing first Chern class and is equipped with a 
(O,l)-connection D = d + a, where a is a (0,l)-form with values in the Lie algebra 
of the complexification of the gauge group. From E we wish to construct a bundle 
E over some portion U of conformally compactified, complexified space-time (say an 
affine patch of some real slice), so our twistor bundle (E,D) will be defined on the 
region U C CP 3 swept out in twistor space by the lines X corresponding to points 

x e u. 

— 2 

We do not assume that D = on twistor space, and so our bundle has an 

2 

almost complex structure. Nevertheless, D necessarily does vanish when restricted 
to any X = CP 1 , because D is a (0, 2)-form. Any complex bundle on CP 1 that 
is topologically trivial is also holomorphically trivial for small enough a, so E\x is 
holomorphically trivial. Therefore we can find N linearly independent global holo- 
morphic sections of E\% that by Liouville's theorem are unique up to a constant 
GL(iV) transformation. The space T(X, E\x) of such sections is a copy of associ- 
ated to the point x in space-time, and so forms the fibre of a space-time Yang-Mills 
bundle E. 
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More explicitly, whenever E\x is holomorphically trivial we can find a smooth 
gauge transformation H(X, A) such that iJ _1 (<9 + a)\xH = d\x, or equivalently 

a\ x = -dH H' 1 . (2.1) 

Here, A is a coordinate on the Riemann sphere X and the (^-operator in (2.1) acts 
in the A direction. As the line X varies in U, the family H (x, A, A) of such gauge 
transformations varies smoothly with x 6 U. The holomorphic frame H := H(x, A, A) 
on X will play an important role in what follows, and exists even for an almost 
complex bundle E. It is unique up to 

H(x,\,\)^H(x,\,\)g(x), (2.2) 

where g depends only on x because it must be globally holomorphic in A to pre- 
serve (2.1), and Liouville's theorem states that any globally holomorphic function on 
a Riemann sphere is constant. 

Although we have constructed a space-time bundle from an arbitrary almost 
complex bundle E on twistor space (subject only to Ci(E) = and the smallness 
of a), the D-operator on twistor space (depending on six real dimensions) cannot 
in general be encoded in a connection V on four-dimensional space-time without 
imposing further conditions. Nevertheless, we will see that the concept of parallel 
propagation along null rays survives, so even in the most general case we can define 
the holonomy around a null polygon, and hence the Wilson loop required for the 
amplitude discussion. In certain circumstances, say for U a Euclidean real slice [21], 
the condition for D to give rise to a generic off-shell connection on space-time can 

— 2 

be characterised by the vanishing of some components of D . (The vanishing of all 
components of D would imply that V is self-dual.) When the twistor D-operator 
does determine a connection on space-time, the Wilson loop defined below agrees 
with that for the space-time connection. We give an explicit demonstration of how 
this works in section 6. 

2.2 Parallel propagation and Wilson loops 

The standard Wilson loop for a curve C in space-time is computed by first identifying 
the holonomy matrix obtained by parallel propagation of the fibre E p of E at some 
basepoint peC around C and back to p again; one then takes its trace. The parallel 
propagation can be represented as the path-ordered exponential integral 

W[A;C] = TrPexp \- J Aj (2.3) 

of the connection A around C, and is a gauge invariant functional of both A and C. 

In the context of scattering amplitudes, we are particularly interested in parallel 
transport around a piecewise null polygon with vertices at points x\ for % — 1, . . . n. 
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The null ray segments correspond to the massless momenta in the scattering process 
via the formula 

Xi - x i+1 = AjAj, (2.4) 



so that the direction of the null ray is determined by the spinor Aj. In twistor space, 
this null polygon corresponds to a nodal curve whose components are holomorphic 
lines Xj (representing the vertices Xi) and whose nodes occur at locations Zi repre- 
senting the null rays (see figure 1). 

Parallel propagation along a null geodesic has a particularly straightforward 
interpretation in twistor space, because the null ray from Xi to Xj+i corresponds to a 
single point Zi in twistor space where the lines Xi and Xi + \ intersect. The parallel 
propagator 



from Xi to Xi + \ tells us how to compare the fibre E x% with E Xi+1 . As we explained 
above, these space-time fibres are the spaces of global sections of the twistor bundle 
E restricted to X; and X i+1 , respectively. But since X, D X i+1 = Zi, these sections 
can be compared directly at Zi. Thus parallel propagation along a null geodesic 
identifies the fibres of E along that geodesic with the fibre Ez i of the twistor bundle. 
In formula?, this is achieved by setting 



because H (x, A) is a holomorphic frame for E restricted to X, so (2.6) is a gauge 
transformation that compares the frames of the twistor bundle on the two Riemann 
spheres X, and X*+i at their intersection point Zi = (Aj, /Uj). The off-shell equivalence 
of these two forms of the parallel propagator will be proved in section 6. 

Continuing all around the null polygon and taking the trace, we find that the 
Wilson loop is represented on twistor space by 



where P Y\ denotes an ordered product, with % increasing towards the left. The first 
line of (2.7) follows directly from concatenating the parallel propagators (2.6) to 
give the holonomy. The second line of (2.7) re-partitions the product in terms of 
the holomorphic frames that live on each line. Recalling that H(xi,X) was defined 
in (2.1) to satisfy 




(2.5) 



U(x i+ i, xf, Xi) = H(x i+1 , Xi) H(xi, Xi) , 



(2.6) 




(2.7) 



(d + a)\ Xi H = 



(2.8) 
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we see that H(xi, X)H(xi, A.j_i) _1 is the solution to (2.8) that obeys the boundary 
condition that it is the identity matrix at Z — Z^\. 

To use (2.7) in a correlation function, we will need to express the holomorphic 
frame H(x{, X)H(xi, Aj-i)" 1 in terms of the twistor field a that will appear in the 
action. This may be done as follows. Let be the inverse d-bar operator on the 
line Xj, with the boundary condition that it vanishes at the point Zi-\. Explicitly, 
if we parametrise the line Xj by 

Z(s)i = sZ { _ x + Z t , s E C. (2.9) 

then for any (0, l)-form uj on Xj, then the function d~ l u) may be defined by 

{8^u){s) = [ G(s,s')Au(s'), (2.10) 

where the Greens function 



x, 



(and hence d i 1 tu) vanishes at Zj_i (where s = oo). Solving (2.8) perturbatively 
shows that H(xi, X)H(xi, Aj-i)" 1 is given by the infinite series 



oo 1 

U d+a 



(2.12) 



where again we have to take care that the terms in (d i 1 a) 1 are appropriately ordered. 
Thus we obtain the perturbative expansion 



n I oo 



W[C n ] =TrPjI [^(d-HZisW* (2.13) 



i=l \i x =o 



for the null polygonal Wilson loop in terms of the twistor field a. Exactly as on 
space-time, computing the correlation function of this operator using the twistor 
action for M = 4 SYM gives the n-particle MHV amplitudes, as we shall see below. 

2.3 A supersymmetric Wilson Loop 

J\f — 4 SYM naturally lives on supertwistor space CP 3 ' 4 , rather than on CP 3 . Thus, 
instead of a connection (0,l)-form a(Z) we really have an M = 4 superfield 

■y ^abcd ^abcd 

A{Z,x) = a(Z) + Xal a (Z) + -XaXb<t> ab {Z) + -^-XaXbXcld + -^-XaXbXcXd g(Z) 

(2.14) 

that may be thought of as a connection (0,l)-form on a bundle over CP 3 ' 4 . The lowest 
component of A is the non-supersymmetric field a, while the Grassmann expansion 
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of A(Z, x) provides the rest of the supermultiplet. As we will see in the next section, 
the twistor action for M = 4 SYM can be written in terms of the superfield A, fully 
off-shell. 

From this perspective, it is unnatural to define the Wilson loop purely in terms 
of the non-supersymmetric field a as in (2.13). Instead, we replace a by the superfield 
A to find an M = 4 supersymmetric Wilson loop on twistor space (that we also call 
W[C n ]): 

n / oo \ 

Wig^TrPn ^(^(WH • (2.15) 
i=l \/ 1= o / 

It is this operator that we use in the rest of the paper. The correlation function of 
this supersymmetric Wilson loop will provide us with the complete planar S-matrix 
of the M = 4 theory, including all N fc MHV amplitudes. We return to the definition 
of this Wilson-loop in space-time in section 6. 

3. The Twistor Action for M = 4 SYM 

M = 4 super Yang-Mills may be defined on twistor space by the action [22] 

S = S 1 + S 2 , (3.1) 
where Si is a holomorphic Chern- Simons action [39, 40] 

S, [A] = [ D 3|4 Z A Tr [A BA + \a*) (3.2) 

for the connection (0,l)-form A on CP 3 ' 4 , and S 2 is a non-local term 

S 2 [A] = g 2 J d 4|8 x lndet V\ x (3.3) 

consisting of the integral of the logarithm of the determinant 1 of the d-bar operator 
T> = 8 + A, restricted to a line X C CP 3 ' 4 . This object is then integrated over a con- 
tour T in the space of lines, corresponding to a real slice of compactified, complexified 
space-time. Thus the action is defined on the space of connection (0,l)-forms A on 
a complex bundle E — > CP 3 ' 4 that is trivial upon restriction to every line in T. See 
[22] for further discussion. 

Taken alone, the field equations of Si state that T> 2 = F 0,2 = 0, so that E 
becomes a holomorphic bundle. The Penrose- Ward transform [42] states that such 

1 This determinant can be understood either via the Quillen construction [41], or, as in twistor- 
string theory [40, 21, 22], as the partition function of a chiral free fermion CFT on X. Technically, 
the 2?-operator in this operator must be taken to act on sections of E\x <£> Cx(— 1) so that it has 
vanishing index and the determinant exists. Physically, the twisting accounts for the spins of the 
fermions on X, and ensures the CFT has no axial anomaly. 
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holomorphic bundles on twistor space are in one-to-one correspondence with self-dual 
Yang-Mills bundles on space-time. The supersymmetric extension provides the rest 
of the multiplet (including the anti self-dual part of the Yang-Mills field) as linearized 
fields coupled to this self-dual background. The role of S2 is to promote the action to 
describe full, rather than anti self-dual, M = 4 SYM (at least perturbatively). This 
can be seen by expanding lndetP|x as the infinite series 

00 1 f 

lndet t>\ = Tr(ln<9| ) + V - / Tr (<TUi <TU 2 ■ ■ • <TU m ) , (3.4) 

where Ai is the field A inserted at some point pi G X. These insertions are connected 
together using the free fermion propagator, or equivalently, the inverse Cauchy- 
Riemann operator d~ l acting on sections of Ox{— 1)- Each term 2 in this expansion is 
an m-particle MHV vertex, while the holomorphic Chern-Simons term contains the 
remaining 3 particle MHV vertex. 

Both the action and the Wilson loop (2.15) are invariant 3 under general complex 
gauge transformations 

V m- g~ l Vg, (3.5) 

where g(Z, Z) is a smooth map from twistor space to the gauge group. Twistor space 
has six (bosonic real) dimensions, so this is a considerably greater gauge freedom 
than in space-time. As explained in [22], this gauge freedom may be exploited to 
relate the twistor action either to the MHV diagram formalism or to the standard 
space-time action for M = 4 SYM. In the next section, we pick an axial gauge 
in which it is straightforward to evaluate the correlation function (1.1), obtaining 
a correspondence with the MHV diagram formalism. On the other hand, one can 
reduce to space-time by imposing the partial gauge fixing that A\x is harmonic with 
respect to an arbitrary Hermitian metric on the Riemann sphere. In this gauge, the 
first two component fields in (2.14) vanish upon restriction to X so that ^4|x is of 
order x 2 an d the expansion (3.4) terminates. S2 then reduces to 

S' 2 = g 2 J d A x Tr ^2G AB G AB + $ ab V Aa V A b + V 6 $ 6c $ cd $ da j (3.6) 

on space-time, where $ is a space-time scalar field, a left Weyl spinor and 
Gab dx AA ' A dx B A , an anti self-dual 2-form corresponding to the twistor fields <fi(Z), 
7(Z) and g(Z), respectively (see [22, 43] for details). Likewise, in this gauge Si may 

2 The term with m = 1 is absent because of the colour trace. 

3 Twistor space has vanishing third Bctti number, so gauge invariance of the holomorphic Chern- 
Simons term is automatic. Because det(T>) is a section of a Quillen determinant line bundle, it is 
well-defined only up to an overall (complex) phase, so lndet (I?) is defined only up to an additive 
piece. This term is annihilated by the fermionic integration in d 4 ' 8 x; see [22]. 
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be reduced to the space-time form 

S[ = J d 4 x Tr \g ab F ab + $ A ' a D AA ,y A a + l -D AA ,® ab D AA '<$> ab + $ a $ A ' a & A , 

(3.7) 

and S[ + 5*2 is the Chalmers-Siegel action for perturbative Af = 4 SYM [34, 33]. 

Because the twistor action is fully off-shell, its associated path integral may be 
used to compute arbitrary correlation functions as well as scattering amplitudes. In 
this paper we focus on the correlator (1.1), but the techniques have much wider 
applicability, as will be explored in [38]. The proposal that the this correlation func- 
tion reproduces the usual scattering amplitudes (divided by the MHV tree amplitude) 
arises if we declare that the action (3.1) and correlator both live in momentum twistor 
space [6, 7]. This is completely analogous to computing the expectation value of a 
polygonal Wilson loop using the standard, space-time M = 4 SYM action and then 
comparing to MHV amplitudes by declaring the computation to have taken place on 
region momentum space (dual space-time). 

3.1 Correlation functions from the twistor action in an axial gauge 

In this section, we perturbatively compute the correlation function (1.1) using the 
twistor action in the axial gauge 

i.4 = 0, (3.8) 



where is a fixed reference twistor. This gauge greatly simplifies the evaluation of 
the correlation function, because the cubic Chern-Simons vertex is eliminated leaving 
just the MHV vertices in S^- However, since these vertices all come with a factor 
of the Yang-Mills coupling g 2 , they are irrelevant when computing the correlator at 
lowest order (g°). Thus, at order g°, the correlator is evaluated 

J[dA] W[C n ] expC-SitA]) (3.9) 

in holomorphic Chern-Simons theory, by pairwise contracting fields in W[C n ] with 
axial gauge propagators A* of the holomorphic Chern-Simons theory. According to 
section 2 we can expand the operator as 

n / oo \ 

W[CJ=TtPn (3-10) 
i=l \k=o J 

where the appropriate cyclic ordering in the matrix products is understood. In the 
correlation function, only terms with an even number of A insertions contribute at 
this order. We will show that the correlator of the m = 2k term in (3.10) computes 
the ratio of the N fe MHV tree amplitude to the MHV tree. 
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At higher orders in the coupling we can also make use of the MHV vertices (3.4) 
in S*2. Each such vertex comes with an integral over the choice of line X C CP 3 ' 4 , or 
equivalently an integral over (dual conformal) space-time. From the point of view of 
scattering amplitudes, these integrals will be the usual loop integrals. The expansion 
of (W[C n ]) in powers of g 2 is thus the loop expansion of the amplitude, with the 
number of loops equal to the number of MHV vertices used in the correlator. Note 
that S' 2 - the space-time equivalent (3.6) of the sum of all MHV vertices - is closely 
related to the chiral Lagrangian insertions used by [37] to obtain the integrand of 
certain loop amplitudes, here in an M = 4 context. 

In this paper, we do not attempt to evaluate the loop integrals themselves - they 
are divergent and require regularisation. However, these integrals can be regulated 
and evaluated directly in twistor space. This has been done in [26, 27] at one loop 
and in [28, 29] for certain two loop integrals, using the Coulomb branch (or AdS) 
regularisation scheme introduced in [30, 44]. Higher-loop MHV diagrams have also 
been considered in [45, 46]. 

In this axial gauge, the Feynman diagrams of the twistor action are MHV di- 
agrams [24, 23], with the CSW reference spinor determined by the fixed reference 
twistor Z* (the reference twistor is usually taken at infinity, where its information 
is that of a spinor). We will see that, remarkably, the space-time MHV diagrams 
generated by the correlator (1.1) are the planar duals of the usual MHV diagrams 
for the corresponding scattering amplitude. 

3.2 The propagator in an axial gauge 

The first ingredient we need is the propagator of the holomorphic Chern-Simons 
theory. In axial gauge this takes the particularly simple form 

A*(Z,Z') = 5 2l4 (Z,Z',Z«) = [ — A5 4 V + c 2 Z' + c 3 Z») (3.11) 

J CF 2 ClC 2 C 3 

where D 2 c = C\dc2 Adc3+ cyclic. This delta function restricts Z and Z 1 to be collinear 
with the reference twistor (in the projective space; see [47, 7] for further discussion 
of such projective 5-functions). To understand this definition, recall [24, 48] that the 
propagator for the A field is the inverse of the <9-operator on CP 3 ' 4 , so may be 
represented by a (0,2)-form A on CP 314 x CP 3 ' 4 that has homogeneity zero in each 
entry and satisfies 

B Z A(Z, Z') = 5 Z ,A(Z, Z') = S^ 4 (Z, Z') , (3.12) 

where the 5-function on the right is a (0, 3)-form supported on the diagonal CP 3 ' 4 C 
CP|! 4 x CP|l 4 . The only possible singularities of <P' 4 (Z, Z' , Z*) occur when any two 
of the three points collide, where being collinear is no longer a restriction. Indeed, a 
direct calculation using the definition (3.11) gives 

dA*(Z, Z') = 5 3 I 4 (Z, Z') + 5 3|4 (Z\ Z*) + 5 3|4 (Z„ Z) . (3.13) 
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The first term is the required singularity of the propagator on the diagonal. The 
remaining two terms are 'spurious' singularities that are a common feature of axial 
gauges. A proper study of this propagator, including its relation to more standard 
momentum space propagators and its use in the MHV diagram formalism for scat- 
tering amplitudes on ordinary twistor space, will be presented in [49]. 

4. Tree amplitudes 

In this section we compute the NMHV and N 2 MHV tree amplitudes from the Wilson 
loop. The answers are given as a sum of products of dual superconformal invariants, 
and agree with the results obtained from the momentum twistor reformulation of 
the MHV rules derived in [1]. There it was shown that the MHV diagram rules 
in momentum twistor space associated just 1 to each vertex, and a dual-conformal 
R-invariant (4.4) to each propagator. Here we will similarly see an R-invariant aris- 
ing for each propagator but this will be associated to the planar dual of the MHV 
diagram. 

4.1 NMHV 

The zeroth order term is just (Tr (Id)). For an SU(iV) gauge group, we divide by N 
to normalize this term to 1. The first nontrivial term to compute is thus 



where each field is inserted at some point on the components Xj and Xj of the curve 
C n . (The integration over Xj and Xj is implicit in the definition of d~ l .) 

Suppose first that % ^ j so that the fields are inserted on different components 
(see figure 2). We may parametrize their locations by 



in terms of local coordinates s and t on the two lines. The operator 0~ 1 |x i was 
defined as the Green function on Xj that vanishes at Z^_i, i.e., at s = oo. Thus the 
Greens functions are simply ds/s and dt/t. Contracting the two fields with a single 
twistor space propagator gives 




(4.1) 



1<J< i<n 



Z(s) = sZ^x + Zi 



and Z(t) = tZj-x + Zj 



(4.2) 




(4.3) 



where the second line follows from the definition (3.11) of the propagator, after a 
rescaling of the integration variables. Without further calculation, we can recog- 




Figure 2: A single propagator connecting two ^4 insertions on lines X, and Xj corresponds 
to the basic dual super confer mal invariant 1, i, j— (The diagram shows only the 
relevant lines of C n . 

nise (4.3) as a dual super confer mal 'R-invariant' in the momentum twistor form 4 

/d 4 t 
— — — <5 4|4 (Z a + tiZ b + t 2 Z c + t 3 Z d + t A Z e ) 



5°l 4 (^ a (6, C ,rf,e) + cycli 



ic 



(4.4) 



(a, 6, c, d) (6, c, d, e) (c, d, e, a) (d, e, a, b) (e, a, 6, c) 
first introduced in [7]. Thus we have immediately 

(Tr (cTU, Br^Aj) = [ * , i- 1, i, j - 1, j] (4.5) 

as the contribution to the correlator from two fields inserted on two different lines, 
as illustrated in figure 2. 

We now sum over the possible locations for the field insertion. The R-invariant 
vanishes when two or more of its arguments are the same because [* , i — 1, — 
is totally antisymmetric and homogeneous of degree zero. Thus the contribution 
vanishes whenever the two fields are inserted either on the same line or on adjacent 
components of C n . Summing over components therefore gives 

<Tr(«rU<rU))= Yl (4-6) 

l<i<j<n 

to order g°, where the inverse d-bar operators here are understood to live on C n {i.e., 
they include the sum over components of C n in their definition). This is the precisely 
the momentum twistor expression 

^Smv = XJ[*,i-l,i,J-l,j]- (4-7) 

i<j 



4 In the notation of [50], [* — — = R*sj. We use the notation [,,,,] to emphasise 
that the R-invariant is a totally antisymmetric function of five supertwistors, homogeneous of degree 
zero in each entry. We also define (a, b, c, d) to be the skew product of the bosonic components of 
four supertwistors Z a , . . . , This notation was introduced in [31]. 
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Figure 3: Correlation function diagrams with two propagators (and no vertices) contribute 
to the N 2 MHV tree amplitude. Only those diagrams whose contractions are compatible 
with the ordering of the trace in the Wilson loop survive in the planar limit. 

for the n-particle NMHV tree amplitude (divided by the MHV tree) found in [1] 
using the MHV diagram formalism. In terms of MHV diagrams for the scattering 
amplitude, each term in this sum comes from a diagram that connects two MHV 
vertices with a single propagator. Notice also that if we choose the reference twistor 
Z* to equal one of the external momentum twistors, say Z n , then (4.6) reduces to 
the standard BCFW form of the NMHV tree [50]. 

4.2 N 2 MHV 

The next non-vanishing term at order g° contains four powers of A. Consider first 
the case where these fields are each inserted on different lines, say 



Z(s) = sZ-i^i + Zi, Z{t) = tZj-! + Zj, Z{u) = uZ k _i + Z k , Z{v) = vZt_i + Z h . 



Although all possible contractions are allowed, only those in which the contracted 
.As are adjacent in the trace ordering survive in the planar limit, the alternating 
contraction being suppressed by a factor of 1/N. Thus, only the diagrams where the 
twistor propagators do not cross survive in the planar limit (see figure 3), although 
we can allow i = I and/or j = k. 

In the generic case i ^ I and j ^ k we find 



where each factor of the integration measure is provided by an inverse <9-operator 
on C n . Once again, the identification as a product of two dual superconformal 
invariants (depending on the axial gauge reference twistor Z#) follows immediately 
from the form (3.11) for the propagator. 



(4.8) 




(4.9) 



* ,i— j — 1, j] [*, k-1, k, 1-1,1} , 
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Finally, there is an exceptional class of diagrams to consider where more than 
one field is inserted on the same twistor line. Suppose first that j = k so that we 
have 



Z(s) = sZ^i+Zi, Z(t{) = ttZj-x+Zj, Z(t 2 ) = t 2 Z k _ 1 +Z k , Z(v) = vZ^+Zi . 

(4.10) 

There is no difference in the form of the propagators A*(Z, Z'), nor in the d^ 1 - 
operators on the lines — and (k — l,k). However, compared to the measure 
in (4.9) we must replace 

dt du dti dti 

> 7 ^ — (4-11) 

because, using the cyclicity of the trace and reading from the right, the <9 _1 operators 
in • • ■ d~ 1 A(s 2 ) B~ 1 A(si) <9 _1 ■ ■ ■ first link the node at Z f (where Sj = oo) to A(s\), 
then link A(si) to A(s 2 ) on the same line component, and then link A(s 2 ) to the 
node at (where s» = 0). Consequently, this term gives a contribution 

-^-^--A*(Z( Sl ),Z(t))A*(Z(u),Z(s 2 )) 

s h t 2 - h v (4.12) 

= [*,k-l,k,i-l,i] 

that differs from (4.9) only by replacing Z k -i by the intersection of the line (k—1, k) 
with the plane ( * , i — 1 , i) : 

z k _ x -> z k ~T 1 = (k-i,k)n(*,i-i,i) 

= (*,i-l,i,k-l)Z k - (*,i—l,i,k)Z k -i . 

This follows from a simple change of variables t 2 — > w = t 2 — t\ in the integrand 
of (4.12), with identified as the coefficient of w inside the 5-function. 

Terms in which both propagators end on the same pair of lines, or where one 
of the propagators has both ends on the same line vanish identically, while terms 
in which the A* propagators connect the fields in a different order (so that s\ and 
s 2 are exchanged in the twistor propagators of (4.12), but not in the measure) are 
suppressed in the planar limit. Therefore, summing over the possible lines on which 
the propagators end gives 

(Tr (<rU<rU<9-U<rU)} = K^^j-M [*,jPi,jm-i,z] (4.i4) 

l<i<j<k<l<n+i 

to order g°, where the hatted variables are defined as 

. (i-lJ) n (*,k-l,k) if / =i , 
i-l= { V ; V ; (4.15) 

i — 1 otherwise, 



- 15 - 




Figure 4: The Feynman diagrams of the correlation function (shown in black) are dual 
to the MHV diagrams for the corresponding amplitude (shown in grey). 

and 

, — ((k-l,k)n(*,i-l,i) rfk = j 

k-1 = \ v ' ' y ' ; J (4.16) 

I k — 1 otherwise. 

This is precisely the momentum twistor form of the N 2 MHV tree amplitude, divided 
by the MHV tree, as found in [1]. From the perspective of MHV diagrams for the 
scattering amplitude, term by term, each summand in (4.14) comes from a diagram 
consisting of two propagators joining three MHV vertices. The shifts (4.15)-(4.16) of 
the momentum twistors occur when two propagators border the same region, so that 
they are adjacent in the cyclic ordering of the middle MHV vertex (see figure 4). 

4.3 Dual MHV diagrams and N fc MHV amplitudes 

It is revealing to draw MHV diagrams corresponding to the twistor configurations for 
the correlation function as found above. Lines in twistor space correspond to points 
in space-time, so we can represent a twistor propagator whose ends are integrated 5 
over the lines i) and by a space-time propagator joining the two points 

Xi and Xj. As shown in figure 4, the resulting diagrams for the correlation function 
are the planar duals of the usual MHV diagrams [24] for the corresponding NMHV 
and N 2 MHV tree amplitudes. 

Dualization leaves the number of propagators unchanged, so at order g°, the 
N fe MHV tree amplitude will come from the expansion of (W[C n ]) to order 2k in 
the field. These fields are contracted using k propagators that, in the planar limit, 
do not cross each other. The propagators therefore divide the polygon into k + 1 
regions, so the dual MHV diagram is a planar diagram with k + 1 MHV vertices 
and k propagators. As described in [1], such MHV diagrams yield a product of k 
R-invariants in momentum twistor space. The boundary terms where two or more 
Wilson loop propagators end on the same edge correspond precisely to the boundary 
terms for the MHV diagrams in momentum twistor space, where two propagators 

5 Given two arbitrary lines in twistor space, there is a unique line that passes through a generic 
point Z. t and intersects both given lines. Since P*(Z,Z') requires Z, Z' and Z* to be collinear, 
it follows that for a generic choice of Z*, only one twistor propagator actually contributes to the 
integrals. 
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are adjacent to the same region. The shifts of the arguments of the R-invariants 
in [1] follow in exactly the same way as for the N 2 MHV term above, and generalise 
to arbitrary N fc MHV amplitudes. It is therefore easy to see that our correspondence 
extends to arbitrary N fc MHV amplitudes. 

It is also easy to see that, had we defined our twistor space Wilson loop using the 
non-supersymmetric connection D = d + a that corresponds to the standard space- 
time Wilson loop, then all the correlation function diagrams discussed in this section 
would vanish. This is because the kinetic term of the holomorphic Chern-Simons 
theory on CP 3 ' 4 has a component expansion 



so the component propagator connects an a field to a g field, rather than connecting 
two a fields. To zeroth order in the Yang-Mills coupling, there would be no way 
to contract any field insertions. With this non-supersymmetric operator, the only 
non-vanishing term would have been the trivial one, representing the MHV tree. 

5. Loop amplitudes 

At higher order in the Yang-Mills coupling, we may connect the field insertions in 
W[C n ] via any of the MHV vertices in the action (3.3). As explained in section 3, 
these vertices are supported on an arbitrary line X C CP 3 ' 4 , unrelated to the compo- 
nents of C n . The choice of X is integrated out over a contour V corresponding to a 
real slice of space-time. These (dual) space-time integrals, coming from the non-local 
term in the twistor action, are the loop integrals of the original amplitude. 

We can rewrite the integral over space-time in a way that emphasises its super- 
conformal invariance. If X is the line (A,B), we have that 



where Xa,b define the location of A and B on (A, B). The rest of the integrand is 
independent of these locations. Hence, if we modify the definition of the contour 
T to include a factor of the anti-diagonal 6 S 2 C CP^ A x CP^ s for each x, then the 
Xa, Xb integrals simply yield 1 and the dependence on the specific reference points 
is removed. This choice of contour is particularly natural in twistor space. Although 
the new contour looks like 1R 4 x S 2 locally, globally the 5' 2 s can be chosen to fibre over 
the space-time base in such way that the whole contour is the anti-diagonal CP 3 C 
CP^ x CP^ (though the integrand requires regularisation on this contour). The 

6 That is, we pick a notion of complex conjugation A — > A and set Xb = Xa- The A-dcpcndcnt 
measure (A^dA^) A (A^dAs) / (XaXb) 2 reduces to the standard Kahler form on this S 2 . 




(4.17) 



B 3 ^Za A D 3 ^Z b = d 4 ' 8 



x A 



(X A dX A ) A (X B dX B ) 
(XaXb) 2 



(5.1) 
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implied real slice of compactified, complexified space-time has topology 5* 4 and may 
be thought of as a rotation of the Euclidean real slice towards Lorentzian signature 
space-time, in accordance with the Feynman ie-prescription. (See [26, 27] for further 
details.) 

5.1 MHV 

The action contains a two-point MHV vertex. This plays no role in the usual MHV 
diagram formalism (even as a vertex) because momentum conservation makes its 
numerator factor vanish. However, the momentum conserving 5-function is absent 
in the correlation function calculation. Consequently, the first case to consider is 

(Tr(<rU<rM)) (5.2) 

which gave the NMHV tree at order g°. At order g 2 , the two fields are each connected 
to the intermediate line (A, B) - the MHV vertex - using A* propagators. If the 
fields are inserted on C n at 

Z{t) = sZ^ x + Zi and Z(t) = tZj_ x + Z i (5.3) 
and the propagators meet the auxiliary line (A, B) at 

Z(ui) = Z A + UiZ B and Z{u 2 ) = Z A + u 2 Z B , (5.4) 
then the contribution to the correlator is 

/ D 3 I%AD 3 I 4 Z B / . dUldU \ 2 A*(Z(s),Z( Ul ))A*(Z(t),Z(u 2 )). (5.5) 

The measure 

dtti du 2 



(5.6) 

(mi - u 2 y 

comes from the d~ l operators in the expansion (3.4) of the MHV vertices in the 
action. It may be understood as a two-term current correlator in the fermionic CFT 
that generates the MHV vertices. 

It is again straightforward to identify the integrated propagators with dual su- 
perconformal invariants, with the result 

J B 3 ^Z A A D 3 I 4 Z B [ * , i - 1, i, A, B'] [ * , j - 1, j, A, B"\ (5.7) 

as shown in figure 5. Just as in the boundary cases of the N 2 MHV tree, the mea- 
sure (5.6) mixes together the twistors Z A and Zb in the arguments of the 5-functions. 
We can again choose a change of variables so that B is shifted as 

B' — (A, B) r\(*,j — l,j) and B" = (A,B) n — . (5.8) 
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Figure 5: A particular configuration contributing to the 1-loop MHV amplitude. 




Figure 6: The space-time diagram for the correlation function is dual to the MHV diagram 
of the scattering amplitude. 

One can check that the two A fields must be inserted on different components of C n 
for a non- vanishing contribution (although adjacent components do now contribute), 
so summing over the possible insertion points gives 

<Tr(9-U<rU))= / D^Z A AD^Z B V [ * , z-1, i, A, B'\ [ * ,j-l,j, A, B"] (5.9) 
at order g 2 . 

The integrand in (5.9) agrees identically with the integrand of the momentum 
twistor form of the n-particle 1-loop MHV amplitude as presented in [1], which was 
proved to be equivalent to the expression for the 1-loop MHV amplitude obtained 
from the standard MHV rules in momentum space [24, 51, 52]. We thus have 

<Tr («rUcrU)) 0(g2) = M« v . (5.10) 

Once again, the summands in (5.9) correspond term by term to 1-loop MHV dia- 
grams for the scattering amplitude. The arbitrary twistor line X corresponds to an 
arbitrary space-time point x, so the space-time diagram corresponding to figure 5 
has a propagator connecting x^ to x and a further propagator connecting x to Xj. As 
shown in figure 6, this diagram is nothing but the planar dual of an MHV diagram 
for the scattering amplitude. 
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Figure 7: A contribution to the 1-loop NMHV amplitude coming from a cubic term in 
the field expansion of the Wilson loop. 

5.2 NMHV 

As a further example, we consider the integrand of the planar 1-loop NMHV ampli- 
tude. Somewhat surprisingly, this amplitude receives contributions from two separate 
terms in the correlator; specifically, we shall show that M^ HV comes from the order 
g 2 contribution to 

(Tr (<rUrU<rU)) + (Tr (9'UrUrUrU)}. (5.11) 

The first term has three field insertions. When these are each on separate lines we 
obtain 

/ D^Z A A D^Z B d^dna ^ d^ 

Jr i u i ~ u 2 ){u 2 - u 3 )(u 3 - Mi) £A Sa 

where 

Z(si) = sxZi-i + Z h Z(s 2 ) = s 2 Z j _ 1 + Zj, Z{s 3 ) = s 3 Z k ^ 1 + Z k (5.13) 

are the locations of the external fields, connected by propagators to points Z{u a ) = 
ZA + u a Z B , for a = 1, 2, 3 respectively. It is clear from the integrand that this term is 
a product of three dual superconformal invariants, each depending on the reference 
twistor and a pair of external twistors, together with some combination of the 
auxiliary twistors Za and Zg. Performing a change of variables and identifying the 
R invariants, we find the contribution 

J D 3 \ A Z A A D^Z b [ * , i - 1, i, A, B'} [ * j - 1, j, A, B"} [ * , k - 1, k, A, B'") (5. 14) 

shown in figure 7. The shifts in loop variables again arise from changing variables in 
the integrand of (5.12) and are defined as 

B' = (AB) n (*, fc-1, k) , B" = (AB) n (*, z-l,i) , B'" = (AB) n (*, ]) . 

(5.15) 
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As at MHV, the corresponding diagram where two (or more) propagators end on the 
same component of C n can be shown to vanish, so in total we find 

/d 3| %AD 3 I% [*,i-i,i,A,B'][*j-l,j,A,B"}[*,k-l,k,A,B'"} 

JV l<i<j<k<n 

(5.16) 

to order g 2 . 

Diagrams that contribute to (Tr (d^Ad^Ad^Ad^A)) at order g 2 have the 
same form as the 1-loop MHV diagram, but with an extra twistor space propagator 
connecting the remaining two fields directly [i.e., without connecting to (A, B)). If 
the external fields are inserted on distinct C n components 

Z(s 1 ) = s 1 Z i _ 1 +Z i Z(s 2 ) = s 2 Z j _i+Z j Z(s 3 ) = s 3 Z k _x+Z k Z(s A ) = s^Z^+Zi 

(5.17) 

as for the N 2 MHV tree, then (suppressing the integration over the twistors associated 
to the loop) these terms give a contribution 



/ 



TT — / Mld \ ^(Z( Sl ),Z( Ul )) A*(Z(s 2 ),Z(u 2 )) A*(Z(s 3 ),Z(s 4 )) 



= [*,i-l,i,A,B] [*,j-l,j,A,B] [*,k-l,k, 1-1,1], 
where 

B = (AB)n(*,j-l,j) and B = (AB) D (*, i) (5.19) 

are the shifted loop variables. 

If either j = k or I = i (or both) so that two of the fields are inserted on the same 
component of C n , the corresponding B~ l \c n propagators connect the fields directly 
rather than via a node, and consequently one must replace 

ds 3 ds 3 dsi dsi . . 

► 7 r or > 7 V 5 - 20 

S3 (S3-S2) Si (s 4 -si) 

in the measure of (5.18). These replacements lead to a corresponding shift in the 
external twistors Zj-\ = Z k -i or Z\_\ = given below. Once again adding up 

all the possible configurations of A insertions one finds 



Tr (VUd-U<rU<rU| c J 

= |d 3 I%AD 3 I% ^[*,tt,i,A,fy[*J-lJ t A,B][*,jPl,k,l-l,, 



(5.21) 
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Figure 8: The two classes of MHV diagram contributing to the planar 1-loop NMHV 
amplitude are dual to the diagrams of the two terms in (5.11) in the expansion of the 
Wilson loop. 



at order g 2 , where the triple summation is over the range l<i<j<k<l<i + n, 
and where the hatted variables are defined as 

. (i — l.i) fl ( * , k— 1, k) ifl = i + n 
i-l = ->\ ' v ' (5.22) 



and 



i — 1 otherwise, 



fc _ <*-i,*)n(.,A,B) ■« = , (523) 

k — 1 otherwise, 



to account for the exceptional classes of diagram. Likewise, the shifted loop variables 
are here defined as 

B = (AB)n(*J-lJ) and B — (AB) fl (*,i — l,i) . (5.24) 

As promised, combining (5.16) & (5.21) gives the momentum twistor 1-loop NMHV 
amplitude, in the form presented in [1]. There, it was proved analytically that (5.16) 
+ (5.21) agrees with the 1-loop NMHV integrand obtained from MHV rules in mo- 
mentum space. That these expression in turn agree with conventional Feynman 
diagrams of the scattering amplitude was proved for any 1-loop amplitude in [53]. 
We have thus proved 

<Tr («rU<rU<rU)} 0(g2) + <Tr («rU<rU<rU<rU)) 0(g2) = m« hv . (5.25) 

Again, the space-time diagram for the correlation function is dual to the MHV dia- 
gram for the 1-loop NMHV scattering amplitude, as shown in figure 8. 

As at order g°, notice that using the M = connection d + a to define the Wilson 
loop rather than the J\f = 4 connection 8 + A would result in the vanishing of these 
contributions to the correlator. For the diagram shown in figure 7, this is because 
the (a, g) propagator would have to link the three external a fields to three g fields 
on the line X, but being an MHV vertex, X only supports a most two g fields. The 
contribution considered in (5.21) would also vanish for the same reason as at tree 
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level. However, the diagram in figure 5 for the 1-loop MHV amplitude remains non- 
vanishing, because the two propagators perfectly match the two g fields supported 
on X. 



5.3 The two loop MHV integrand 

The final example we shall consider in this paper is the integrand of the planar 2-loop, 
n particle MHV amplitude. This arises from the same two terms 

(Tr (<rU<rU<rU)) + (Tr (<9-M<rU<rU<9-U)) (5.26) 

as the 1-loop NMHV amplitude, but at order g 4 . These diagrams thus involve two 
auxiliary twistor lines, {A,B) and (C,D). 

Contributions from (Tr («9~U <9~U <9~M))) come from using twistor propagators 
to connect two of the field insertions to one of the auxiliary lines and the remaining 
field to the remaining auxiliary line, before finally connecting (A, B) to (C, D) with a 
further twistor propagator. The space-time diagram in shown in figure 9 and is again 
dual to one of the two classes of MHV diagram. Assuming the field insertions to be 
on the separate C n components and (k—1, k) we find a contribution 

f TT ds a yr du b -A- dv c 

(A,B) J (C,D) J S a K - U b+l ) Al (V c - V c+1 ) 

A*{Z{ Sl ),Z( Ul )) A*{Z{s 2 ),Z{u 2 )) A*(Z(u 3 ),Z( Vl )) A*{Z(v 2 ),Z{s 3 )) 

(5.27) 

where 

Z(si) = syZ^ + Zi Z(s 2 ) = s 2 Z j _ 1 + Zj Z(s 3 ) = s 3 Z k _x + Z k (5.28) 

parametrize the locations of the external fields and for b = 1, 2, 3 and c = 1, 2 we let 

Z{u b ) = Z A + u b Z B and Z{v c ) = Z c + v c Z D (5.29) 

parametrize the locations that the propagators are attached to the auxiliary lines. 
We have used the shorthand L AB \ to indicate the integral J r D 3 ' 4 ^ A D 3 I 4 Zb. In 
the case that two external fields are inserted on the same component of C n , say 
= (k—1, k), the different structure of the <9 _1 |(7 n operator with respect to the 
nodes leads to the replacement 

> 7 7 ( 5 - 30 ) 

s\ [si - s 3 ) 

as in earlier examples. This again has the effect of shifting in the arguments 
of the dual superconformal invariants. Summing over all possible insertions, we find 
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+ 




Figure 9: Diagrams contributing to the planar 2-loop MHV amplitude. 

that 

(Tr («rU<rU<rU)) = 

T f I [*,i^l,i,A,B'][*J-lJ,A,B"][*,A,B"\C,D"][*,k^l,k,C,D'\, 

J(A,B) J(C,D) 

(5.31) 

where the sum is over the range \ <i<j<k<i + n and the shifted loop and 
external twistors are 



B' = (A, B) n ( * , C, D) 

B" = (A,B) n 

B'" =(A,B)n(*,j-l,j) 

while the external shifts are 



D' = (C,D)D(*,A,B) 

D" = (C,D)n(*,k-l,k) (5.32) 



%-\ 



k-1 



(i-l,i)n(*,C,D) i£k = i 
i — 1 otherwise 

(k-l,k)n(*,A,B) Hj = k 
k — 1 otherwise. 



(5.33) 



The remaining term (Tr (B~ 1 Ad~ 1 Ad~ 1 Ad~ 1 A)) is essentially just two copies of 
the 1-loop MHV expression, in that two of the fields are attached to the line (A, B) 
while, quite independently, the other two are attached to (C, D) (see figure 9 for the 
space-time MHV diagram). Therefore, in the case that all the external fields are 
inserted on different components of C n , we obtain a contribution 

/ [*,i-l,i,A,B'][*j-l,j,A,B"]x [ [ * , k—1, k, C, D'][*, 1, C,D"] (5.34) 

J(A,B) J (CD) 



that is just the product of two independent copies of a 1-loop MHV contribution. 
The two terms mix when either j = k or I = i, with the familiar effect that Z^_\ 
and/or become shifted. Summing over all possible insertion points leads to the 
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order g 4 result 

(Tr(<9-U<9-U<rU<9-U)) = 

W / [*,i^l ) i ) A ) B'][* )3 -l )3) A,B"][*,k^l 1 k,C,D'][*,l-l,l,C,D"} 

J{A,B) J(C,D) 

(5.35) 

with the shifted loop variables 

B> =(A,B)n{*,j-l,j) D'= {C,D)n{*, 1-1,1) 

B" =(A,B)n(*,i-l,i) D" = (C,D)n(*,k — l,k) 

and external shifts 

'(i-l,i)n(*,C,D) i£l=i 



i-1 



k-l 



i — 1 otherwise 

(5.37) 

{k-l,k)n{*,A,B) ifk = j 
k — 1 otherwise. 



The summation range in this second term isl<z<j<A;</<z, again understood 
mod n. Combining equations (5.31) & (5.35) gives the complete integrand of the 
planar 2-loop n particle MHV amplitude in the form found from the momentum 
twistor MHV rules [1]. 



6. The Correlation Function in Space-Time 

Since our holomorphic Wilson-loop operator is essentially the twistor space tran- 
scription of the space-time Wilson loop, it is clear that if the space-time Wilson 
loop were defined fully supersymmetrically, its correlation function would also give 
all N fc MHV amplitudes. The only question, therefore, is what the supersymmetric 
extension of the space-time Wilson loop should actually be. Clearly, it must com- 
pute holonomy of some superconnection around a null polygon in superspace, with 
vertices at (xj, 9i). In this section we see that the required superconnection is deter- 
mined by a supersymmetric version of the Penrose- Ward transform. This can also 
be used to identify the space-time and twistor space definition of the Wilson loops. 
Note that the twistor action is most closely related to the space-time M = 4 SYM 
action of Chalmers and Siegel [33, 34], so we expect that the resulting correlation 
function should be calculated using this action. We have yet to check whether the 
space-time correlator reproduces the scattering amplitudes. 

6.1 Construction of a space-time superconnection from the twistor data 

To construct a space-time superconnection from the twistor field A, we start from 
the M = 4 generalisation of (2.8), namely 

(B + A)\ x R(x,9,X) = (6.1) 
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that follows because again the supertwistor bundle E ->■ CP 3 ' 4 is trivial when re- 
stricted to any CP 1 , so we can always find a holomorphic frame H on X and, for 
small enough A, we can find a family H(x,0,A) of such frames simultaneously for 
every line X in some region of CP 3 ' 4 . 

In this equation, the restriction of P to X is implemented by imposing the 
supersymmetric incidence relations 

» A ' = -ix AA ' \ A , Xa = e A X A . (6.2) 

The remaining variables (x, 6, A) are coordinates on the spin bundle S and geomet- 
rically, the incidence relations define a projection p : S — > U of the spin bundle to 
a region U in twistor space (corresponding to a space-time region U). As the line 
X C CP 3 ' 4 varies, restricting T> to X really means pulling back the twistor bundle 
(E, T>) to S to obtain p*(E, T>). The solution H to (6.1) gives a frame of p*E that is 
holomorphic over the CP 1 fibres. 

The operators X A d/dx AA \ X A d/d8 A , d/dX annihilate functions on S that are the 
pullback of holomorphic functions on twistor space, since these depend on (x, 9) only 
via the incidence relations (6.2). We can use the pullback of T> to extend these 
operators to the covariant operators 



and 



V- x :=-t+A> (6.4) 

that act on sections and frames of p*E. In these operators, A a 1 and Aq are the 
horizontal and vertical components of p* A with respect to the fibration p : S — > U. 
That is, 

d d 
Aa ' = xA d^ J P* A and Ao = -^jp*A. (6.5) 

In terms of these operatosr, (6.1) is the equation T>jfl = on S. 

Although {X A d/dx AA \ X A d/d8 A , d/dX} commute and form an integrable distri- 
bution, the covariant operators do not obey 



0, pD> A /,B a ] = 0, [D a ,D b ] + = 0, p,V- x ] = (6.6) 



unless the twistor P-operator itself obeys T> 2 = 0. This condition would imply that 
E was a holomorphic bundle of twistor space, corresponding to a self-dual connection 
on space-time as in the standard Penrose- Ward construction. Clearly, we do not wish 
to impose this in the off-shell context of a correlation function. However, without 
imposing any conditions, T> does not give rise to a (super) connection on space-time 
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at all, because it depends smoothly on six real variables rather than four. The 
appropriate condition that gives an off-shell space-time connection is 

p,V- x ]=0. (6.7) 

To see this, note that when (6.7) holds we can act on (6.1) with D to give 

V- X (BR) = . (6.8) 

This in turn implies 

9(H- 1 DH) = 0, (6.9) 

where d is the d-bar operator acting on the As. Thus, although H(x, 8, A) only 
depends smoothly on A, the combination H _1 DH is in fact holomorphic globally on 
the A-Riemann sphere. Since it is clearly homogeneous of degree +1 (because D is), 
an extension of Liouville's theorem says that H _1 OH is in fact linear in A, so that 
we have 

IT^OH = {\ A A AA ,{x, 6) , X A T a A {x, 6)) (6.10) 
where (Aaa',T a ) depend only on (x, 6). Thus we have constructed a superconnection 

Daa> = - Aw , Dl = ^X-r a A (6.H) 

on the bundle E' over chiral super-space-time. 

Rearranging (6.10) and using the definition (6.3) we see that H(x,8, A) satisfies 

\ A D AA ,E = -A A 'ft and X A D a A E = 0. (6.12) 

These equations imply that the space-time superconnection satisfies the integrability 
conditions 

\P(Ai Db)b>] = , [Dfr,D%\ + = 0. (6.13) 

While these conditions constrain the form of the component expansion of the super- 
connection, they do not imply the field equations, so the superconnection is off-shell. 

Equation (6.12) also allows us to make contact with the parallel propagator along 
null geodesies and hence Wilson loops in Minkowski signature. If U is Minkowski 
space-time, then fibres of S over twistor space are the null geodesies. Therefore any 
form on S that is pulled back from a form on twistor space will vanish when con- 
tracted into the vector field X A X A 'd A A'- In particular X a 'Aa' = 0, so that contracting 
the first equation in (6.12) with X A ' we find 

X a X a 'Daa'H = , X A D a A E = . (6.14) 

These equations uniquely fix H to be the parallel propagator along (super) null 
geodesies. Explicitly, if 7 : [0, 1] — > U is a curve in chiral space-time with end points 
7(0) = (xi,9i) and 7(1) = (x i+ i,9 i+ i) and tangent vector 

4l)= A ^'^ + A ' 4 ' te SF (615 » 
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X j 

Figure 10: The integrablity conditions (6.13) imply that the curvature of the space-time 
superconnection vanishes above any null geodesic in the non-supersymmetric space-time. 
The Wilson loop is thus independent of the choice of lift of this null geodesic, provided its 
has the same endpoints. 

(t G [0, 1] being a parameter along the interval), then integrating (6.14) gives 
K- 1 {x i + 1 ,d iH] X i )lI(x i ,e i ',\) = Pexp (-J j*(A AA ,dx AA ' + r/d<)J (6.16) 

as a supersymmetrization of (2.6). Repeating for each leg of the Wilson loop will lead 
to the space-time formulation of the supersymmetric Wilson loop denned in (2.15) 
for twistor space 7 . We expect that the correlation function of the space-time su- 
persymmetric Wilson loop arising from (6.16) also computes the complete planar 
S-matrix of M = 4 SYM. 

Note that although the relation to scattering amplitudes provides us with a 
preferred tangent vector (6.15), the integrability conditions 

X A X B [D A a , D B b'] = , \ A \ B [D A a , D B b ] + = (6.17) 

of equation (6.13) imply that the curvature of the superconnection vanishes provided 
we remain above any null ray in the bosonic space-time. Therefore, at least at the 
classical level, the parallel propagator (6.16) is independent of the choice of lift of the 

7 When U is a real 4-dimensional slice of complex conformal Minkowski space that admits no null 
tangent vectors, U has 6 real bosonic dimensions in twistor space and the condition (6.7) leads to 
the vanishing of two of the three components of T> 2 . The reconstruction of an off-shell connection 
is then straightforward. However, when U is a real slice of Lorentz signature, U degenerates to 
become a 5-dimcnsional CR manifold, PN, defined by Z ■ Z = 0. On such a CR manifold, T> 2 has 
only one component and so (6.7) implies T> 2 = 0. Thus, if T> defines a connection on space-time, it 
is necessarily self-dual. In order to obtain the above identification of twistor space with space-time 
Wilson loops for a completely off-shell connection, in this situation we must use a limiting argument 
from slices with no null tangent vectors analogous to standard arguments for obtaining correlation 
functions for operators inserted at light-like separated points obtained as limits of those inserted at 
space-like separated points. 
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null geodesic XiXid/dx to superspace (see figure 10). Said differently, there exists a 
gauge for (A,T) in which the dependence of (6.16) on the Grassmann coordinates 9 
comes only from the endpoints. 

6.2 Component form of the space-time superconnection 

Finally, as an exercise we briefly indicate how the component expansion of this space- 
time superconnection is related to the component expansion of the twistor M = 4 
superfield. This is again determined by the integrability conditions (6.13) and may 
be obtained explicitly be expanding out all the fermionic coordinates at each stage 
in the above procedure. We calculate the superconnection to order 9 2 in the general 
case, and give the full answer at the end for the abelian case. 
The twistor superfield is 

j ^abcd ^abcd 

A(Z, x) = a(Z) + Xa l a (Z) + -XaXb (f) ab (Z) + -^-XaXbXc Id + -^-XaXbXcXd g(Z) , 

(6.18) 

where the component of (xY is a (0,l)-form on the appropriate region U C CP 3 , 
homogeneous of degree — r. Under the Penrose-Ward transform, 7 a (Z) and 'ya(Z) 
correspond respectively to positive and negative helicity gluinos, <fi ab (Z) to scalars, 
and g(Z) to the anti-self-dual part of the Yang-Mills field. We can similarly expand 
equation (6.1) in the fermionic coordinates 9 and solve term by term. Set 

R(x, 9, A) = H(x, A) + 9 A h a A (x, A) + 9 A 9*h a A b B (x, A) + . . . + 6 8 h 8 (x, A) (6.19) 

and expand out equation (6.1) using this, (6.18) and (6.2) for x a to get the component 
equations 

= (B + a)\ x H(x,X) 

= (d + a) \xh A (x, A) + X A 7 a (x, X)H(x, A) 

= (8 + a) \ x h%{x, A) + X A 7 a (x, X)h b B (x, A) + X A X B ab (x, X)H(x, A) (6 ' 20) 



where the last equation is understood to have the symmetries of 6^6^ imposed on 
it. 

The leading order part is the non-supersymmetric form (2.8) that fixes the holo- 
morphic frame H(x,X), while the subsequent terms are what one expects for the 
Penrose transform for fields that are background coupled to the connection arising 
from the first term. We can solve these component equations hierarchically. The 
first leads to the ordinary space-time connection A(x) = A(x, 9)\g=o, because we use 
the bosonic part of (6.10) at zeroth order in 9 to fix A(x) in terms of H as 

X A A AA (x) = X A H~ 1 d AA H . (6.21) 
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The subleading parts can be simplified by using if as a gauge transformation to 
replace the d + a|X operator simply by d. The second equation of (6.20) becomes 

B{H~ l h a A (x,X)) + X A H-^ a H = 0, (6.22) 

where here and in what follows, d always refers to the d-bar operator on X. The 
related equation 

dip a + if-V(x,A)i? = (6.23) 
has a unique solution ip a that has weight —1 in A. The quantity 

:= \ A D AA ,r (6.24) 

is holomorphic in A because X a Daa' annihilates H~ x ^ a H^ and since ^ has weight 
zero in A, Liouville's theorem states that it is in fact independent of A. Therefore, 
setting 

E~ x h\ = i) a \ A (6.25) 

and using (6.10) we determine the coefficient of 9 in the expansion of the connection 
to be: 

A AA ,(x, 9) = A AA ,{x) + 9 Aa ^ a A ,{x) + 0{9 2 ) , T a A = + 0(9) (6.26) 

Note that we are free to add a function p A (x) that is independent of A to the solution 
of (6.22). Doing so leads to the addition of a term D AA >(p b B 9 B ) to A AA i and p a A to 
T A , representing a supersymmetric gauge transformation. 

We can proceed similarly with the third equation in (6.20). Using the solutions 
of the first two equations in (6.20), we may write the third as 



= diH-'h^ix, A)) + \ A (H~ 1 'y a H) (H' 1 h b B ) + XAbH^^H 
= d{H- x h%{x, A)) - \ A \ B {Br)^ b + XaXbH-^H . 
As with the gluino, we first consider the related equation 

d(H~ l h a A \x, A)) + X A {d4j a )ij b + X A H~ 1 <p ab H = (6.28) 

which has a unique soluton H~ l h A = m A (x, A) of homogeneity —1 in A. We can 
now observe that 

: = \ A m A b (6.29) 

has homogeneity degree zero and is holomorphic on the sphere and so depends only 
on x. We identify it with the scalar field in the multiplet. Setting H~ l h a AB = m ab A XB) 
we now observe that 

d{X B D BA ,m a A b + X A 4j [a V b l) = (6.30) 
so by a now familiar argument 

Q% := X B D BA ,mi + X A ^% (6.31) 
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depends only on x and by contraction with X A we discover that 

= D AA ,<S> ab (6.32) 

and we find from equation (6.10) that 

Aaa'bc = e A{B D C )A^ ab (6.33) 

is the coefficient of 9 2 in the expansion of A(x,9). Therefore, to second order in 9, 
we have the super connect ion 

A AA/ (x, 9) = A AA ,(x) + 9 Aa V A ,(x) + 9 Aa 9 B D BA ,$ ab + 0(9 3 ) 
T a A (x,9) = 9 Ab <S> ab + 0(9 2 ) 

Proceeding in a similar way (see also [54]) will lead to a complete component expan- 
sion of the space-time superconnection in terms of (derivatives of) the component 
fields. This is easily performed in the abelian case where we obtain 



A AA ,{x, 9) = A AAI {x) + V a A ,(x)9 Aa + D BA ,$ ab 9 2B ab + D BA ^ Cc 9 A BCc + D BA ,G CD 9 A Ba 
T\( Xl 9) = <b ab 9 Ab + V^aL + GBce f 



(6.35) 

where we have introduced the notation 

2AB q(AqB) niABCa ^abcdnAnBnC MABCD ^abcdnAnB nC nD 



= 9^6 b > , 9 6ABVa = 6^9*9^ , 9 AABUU = e abcd 9*9 B 9^9 u d (6.36) 
(all symmetric on their spinor indices). 
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